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We present an experimental and theoretical study of the polarized photoluminescence spectrum 
of single semiconductor quantum dots in various charge states. We compare our high resolution 
polarization sensitive spectral measurements with a new many-carrier theoretical model, which 
was developed for this purpose. The model considers both the isotropic and anisotropic exchange 
interactions between all participating electron-hole pairs. With this addition, we calculate both 
the energies and polarizations of all optical transitions between collective, quantum dot confined 
charge carrier states. We succeed in identifying most of the measured spectral lines. In particular, 
the lines resulting from singly-, doubly- and triply- negatively charged excitons and biexcitons. 
We demonstrate that lines emanating from evenly charged states are linearly polarized. Their 
polarization direction does not necessarily coincide with the traditional crystallographic direction. 
It depends on the shells of the single carriers, which participate in the recombination process. 

PACS numbers: 73.21. La, 78.67.Hc 



I. INTRODUCTION 



Quantum dots (QDs) are nano-structures, which con- 
fine electrons and holes in all 3 dimensions. This con- 
finement results in a discrete spectrum of single carrier 
energy levels and spectrally sharp optical transitions be- 
tween them. The photoluminescence (PL) spectrum of 
single self-assembled semiconductor QDs is usually com- 
posed of many discrete spectral lines. The variety of 
lines originates from optical transitions between vari- 
ous many carrier configurations and different QD charge 
statesi££ii£. 

Several experimental techniques are used for identifying a 
given spectral line by associating it with a specific optical 
transition. These techniques include excitation intensity 
dependent PL spectroscopy, which distinguishes between 
single-exciton and multi-exciton transitions 6 and second 
order intensity cross-correlation measurements, which de- 
termines the temporal sequence by which the emission 
occurs in general 6 -, and identifies radiative cascades in 
particular— PL excitation (PLE)^ 9 - as well as electro^- 
and magneto 11 ' 12 PL spectroscopies are used to further 
provide information regarding the QD's charge state dur- 
ing the optical transitions. 

Unfortunately, even when all of these methods are com- 
bined, occasionally, some lines still remain unidentified^. 
Polarization sensitive PL and PLE spectroscopy have also 
been applied to aid in line identification. Most notably, 
the neutral exciton and neutral biexciton lines are split 
into two cross linearly polarized doublet o 14 ' 15 : 16 , while 
singly charged excitonic lines are unpolarized, and dis- 
play large circular polarization memor y 9 : 17 : 18 when ex- 
cited quasi-resonantly. 

In this work we focus our studies on polarization sensi- 
tive PL spectroscopy of single semiconductor quantum 
dots. We carefully measure the polarization of the PL 



spectra under various excitation conditions. Our results 
are then compared with, and analyzed by, a novel theo- 
retical many charge-carriers model. The method used 
for the calculation of the many-carrier states and op- 
tical transitions between them is the full-configuration- 
interaction (FCI) method^ 9 -. The novelty in our model is 
in its inclusion of the electron-hole exchange interaction 
(EHEI)^^^. We show that the model provides a very 
good understanding of the experimental measurements. 
The manuscript is organized as follows: In section HT1 we 
describe the sample and the experimental setup used for 
the polarization sensitive PL spectroscopy. In section [TTT1 
we describe the theoretical model and in section IIVI we 
compare theoretical and experimental results. A short 
summary is given in section fVl 



II. EXPERIMENTAL METHODS 
A. Sample 

The sample was grown by molecular beam epitaxy on 
a [001] oriented GaAs substrate. One layer of strain- 
induced InGaAs QDs was deposited in the center of a 
one wavelength GaAs spacer layer. The height and com- 
position of the QDs were controlled by partially covering 
the InAs QDs by a 30A thick layer of GaAs and by subse- 
quent 30 seconds growth interruption 22 to allow diffusion 
of In(Ga) atoms from (into) the strained islands. The 
growth resulted in In x Gai_ x As QDs whose exact shape, 
lateral size, composition and strain profile are unknown. 
The sample was not rotated during the growth of the 
strained layer, resulting in a variable density of QDs 
across the sample's surface 6 . The estimated density of 
QDs in the areas that were measured is 10 8 cm -2 . 
The optical microcavity was formed by distributed Bragg 
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reflecting (DBR) stacks of 25 and 11 periods of alter- 
nating AlAs/GaAs quarter wavelength layers below and 
above a GaAs spacer layer, respectively, giving a Q-factor 
of ~500. The spacer layer was grown to a width close 
to the wavelength in matter of the light emitted from 
the QDs due to ground state e-h pair recombinations 
(1 A cavity). The microcavity improves the efficiency of 
photon collection, but limits the energy in which photon 
collection is possible. In particular, emission of photons 
with energies smaller than the microcavity mode energy 
is forbidden. Therefore, the density of QDs which emit 
efficiently is roughly two orders of magnitude lower than 
their actual density 2 ^. In order to electrically charge the 
QDs, a p-i-n junction was formed by n-doping the sub- 
strate and the bottom DBR and p-doping the top DBR, 
while leaving the GaAs spacer intrinsic. An extra AlAs 
barrier was grown inside the GaAs spacer between the 
p-type region and the QDs. This barrier prolongs the 
hole's tunneling time into the QDs at forward bias and 
out of them at reverse bias, with respect to the tunneling 
time of the electron. This enables negative charging upon 
forward bias and positive charging upon reverse bias. 
The top electrical contact of the sample was made of 
a semi-transparent layer of indium-tin oxide in order to 
provide optical accessibility. The sample was not pat- 
terned or processed laterally to prevent obscuration of 
the QD emission and its polarization. 



B. Optical characterization 

For the optical measurements we used a diffraction 
limited low temperature confocal optical microscop o 2 ' 24 . 
The sample was mounted on a cold finger of a He-flow 
cryostat, maintaining temperature of about ^20K. A X60 
in-situ microscope objective was used in order to focus 
cw or pulsed laser light at normal incidence on the sam- 
ple surface. The emitted light was collected by the same 
microscope objective. The objective was accurately ma- 
nipulated in three directions using computer-controlled 
motors. The collected light was spatially filtered, dis- 
persed by a 1 meter monochromator and detected by a 
nitrogen-cooled CCD array detector. The system pro- 
vides diffraction-limited spatial resolution, both in the 
excitation and the detection channels and spectral reso- 
lution of about 15 [ieV per one CCD camera pixel. 
The polarization of the emitted light was analyzed by two 
computer controlled liquid crystal variable retarders and 
a linear polarizer in-front of the monochromator. The 
degree of polarization of the emitted light and its po- 
larization state were deduced by six independent mea- 
surements of differently polarized spectra and calcula- 
tion of the Stokes parameters 2 ^. Throughout this work 
we use the symbol H (V) for linear light polarization 
along the [110] ([110]) crystallographic axis of the sam- 
ple. These in-plane orientations are determined by cleav- 
ing. The symbol D=^ (H+V) (D=-^=(H-V)) is used for 
the 45° (-45°) diagonal polarization, while the symbol 



R=-L(H+iV) (L=-^(H-iV)) is used for the right (left) 
hand circular polarization. 

A general state of polarization can be represented as a 
vector inside the Poincare sphere. Figure [1] shows a vec- 
tor in the Poincare sphere and its relation to the shape 
and orientation of the polarization. 
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FIG. 1: (a) The polarization state represented as a vector 
P on the Poincare sphere, (b) The same polarization state 
represented as the loci of points that the electric field of the 
light obtains, during one period, in a plane perpendicular to 
its propagation direction, so. .3, are the experimentally deter- 
mined four Stokes coefficients^. 



In Fig. [5] we show the PL spectra from a single QD as a 
function of the voltage applied to the sample. The QD 
was excited by a cw 1.47 eV Ti-sapphire laser light. The 
current through the device as a function of the voltage is 
also shown. 

The specific structure of our sample is such that at for- 
ward biases (above ~7 volts) the QDs are negatively 
charged as clearly evident by the abrupt step in the emis- 
sion energy. This injection induced charging mechanism 
is similar to that reported earlier^i^. At large reverse bi- 
ases, however, the QD is increasingly positively charged, 
due to vast differences between the tunneling-out rates 
of electrons and holes^^. 

The spectral line identification in Fig. [2] is based on the 
order by which the lines appear and disappear as the 
voltage on the device increases. Information gained from 
excitation intensity dependence PL spectroscopy (not 
shown) and polarization sensitive spectroscopy (see be- 
low) is also used for this purpose. 

In Fig. [3] we present the measured polarization sensitive 
spectra for the bias voltages indicated by horizontal lines 
in Fig. H 

We note here that the spectral shapes of the observed 
negatively charged lines and in particular the fine struc- 
ture components of X~ 2 , XX~ 2 and X~ 3 are similar to 
those observed also in previous work a 4 ' 27 . 
In Fig. OUa), the QD was on average neutral. The neu- 
tral, as well as the singly negatively and singly positively 
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FIG. 2: (a) Measured PL spectra from a single SAQD, as a 
function of the bias on the device. The QD was excited by 
1.47 eV cw laser light. The various spectral lines are labeled 
by X (XX) for single (double) initial e-h pair occupation and 
a superscript which denotes the QD charge state during the 
recombination. The horizontal solid lines mark the voltages 
for which spectra are presented in Fig. [3] (b) The current 
through the device as a function of the bias voltage. 



charged exciton and biexciton spectral lines are observed. 
The corresponding polarization spectra projected on the 
linear H-V and on the linear D-D axes of the poincare 
sphere are shown in Fig.[3jb). The projections are calcu- 
lated by subtracting the two cross-linearly polarized spec- 
tra dividing by their sums. Division by zero is avoided by 
adding a constant equals to the standard deviation of the 
background noise to each spectrum. The projection on 
the R-L axis of the Poincare sphere was zero to within 
our experimental uncertainty (not shown). From these 
two projections, the actual magnitude and direction of 
the linearly polarized lines can be straightforwardly de- 
termined. In Figs.^c) andOJd) we present the PL spec- 
trum and its linear polarization projections, respectively, 
for a bias voltage of 7.15 volts at which the QD was on 
average, negatively charged with two to three electrons. 



III. THEORETICAL MODEL 

The model that we developed is a relatively simple 
many-carrier model which includes the electrostatic in- 
teractions between the QD confined charge carriers. Un- 
like previous, similar models 2,5,19 , which neglected the 
electron- hole exchange interaction (EHEI), our model 
includes it. This interaction is indeed orders of mag- 
nitude smaller than the direct Columbic terms. Spec- 
trally, it is only significant when the fine excitonic struc- 
ture of the spectrum is considered. However, when 
the polarization spectrum is considered, this anisotropic 
interactio n 20 ' 28 ' 29 is by far the leading term. 
Our model is constructed as follows: 
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FIG. 3: (color online) (a) and (c) Measured PL spectrum 
for bias voltage of 0V and 7.15V, respectively. The energy 
is measured from the energy of the X° line, (b) and (d) 
PL polarization spectra for bias voltage of 0V and 7.15V, 
respectively. The black (orange) lines present the polarization 
as projected on the H-V (D-D) axis of the Poincare sphere. 
Vertical dash lines at various spectral lines are drawn to guide 
the eye. 



We first solve a single carrier problem for the electron 
and for the hole in the QD. In this manner, we obtain a 
consistent set of single charge carrier energies and asso- 
ciated wavefunctions. 

We then use this set of energies and wavefunctions in 
order to construct a many-carrier second quantization 
Hamiltonian, which includes the electrostatic interaction 
between the confined carriers. In particular we consider 
the EHEI which is introduced into our model using a 
semi-phenomenological approach. 

The many body Hamiltonian is then diagonalized, thus 
obtaining the collective many carriers energies and wave- 
functions. We then use the dipolc approximation to cal- 
culate the optical transitions between the many carriers 
states for a given light polarization. From these calcula- 
tions we finally construct the polarization sensitive emis- 
sion spectrum, which is used for comparison with the 
experimental measurements. 
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A. The single-carrier problem 

The single-carrier energies and wavefunctions are cal- 
culated using the slowly varying envelope function ap- 
proximation (SVEFA) 30 . We use one (doubly Kramers 
degenerate) band for the electron and one band for the 
hole without band mixing ('One-band SVEFA'). This 
approximation results in two independent Schrodinger 
equations for the envelope functions of the electron and 
that of the hole. The potential of the QD is approximated 
by a finite three dimensional potential well in the form 
of a rectangular slab, with the long (short) side oriented 
along the H (V) direction, and with different dimensions 
and offsets for the two types of carriers. The parameters 
that we used are listed in table Q] We solved the differ- 
ential equations numerically, using the finite differences 
method, thus obtaining the single particle eigen-energies 
and envelope wavefunctions. 



TABLE I: The QD parameters used to calculate the single- 
carrier energies and envelope wavefunctions. 



Parameter 


Value 


Units 


QD shape 


Rectangular slab 




QD size for the electron 






(Length x Width x Height) 


244 x 232 x 34 


A 


Hole/Electron length ratio 


0.72 




Electron effective mass 


0.065^- 


m 


Hole z-direction effective mass 


0.342i 


m 


Hole in-plane effective mass 


0.25 


mo 


Electron potential offset 


-324 


meV 


Hole potential offset 


-108 


meV 


GaAs band gap 


1.51922 


eV 


Band gap of QD material 


1.087 


eV 


Ino.5Gao.5As dielectric constant 


13.822 




E p of Ino.5Gao.5As 


25.52i 


eV 



B. The many-carrier Hamiltonian 

The second quantization many-carrier Hamiltonian for 
the QD containing both electrons and holes is given 
b>Al&: 

H = Ho + H ee + H hh + H eh (1) 

where Hq is the single carrier Hamiltonian, H ee (Hhh) 
is the electron-electron (hole-hole) interaction Hamilto- 
nians, and 

»1>*4 ,32 ,33 

is the electron-hole interaction Hamiltonian. The 
electron creation operator a\ and the hole creation 
operator seperately satisfy the regular Fermionic 
anti-commutation relations. 

The quantities of the form Cf l 1 1 p ^ 3 ^ ni , where P1..4 can 
be either 'e' (for electron) or h' (for hole), and the 
indices 711..4 run over the appropriate states, are the 
Coulomb interaction integrals. 

Computation is only feasible with limited number of 
single carrier states. Therefore, only the first 12 lowest 
energy electron and hole states are usually considered in 
our calculations. 

For the single-carrier wavefunctions which are calculated 
using the SVEFA, the Coulomb integrals can be sepa- 
rated into long-range (inter- unit-cell), and short-range 
(intra-unit-cell) integrals. The long-range integral can 
be expanded into a Taylor series in r\ — r 2 — (R\ — R2), 
where Ri(2) is the lattice vector nearest to the position 
vector ^1(2)- In most cases the zeroth order of that series 
is much larger than other orders and also larger than the 
short range integral^. For the one band SVEFA, the 
zeroth order of the long-range interaction is given by^S: 



C£KU = S P uP*8 P 2, P Js ni ,s n Js n2 ,s n3 I I dVdVaC^C/ (^)^r— rr€l(r2)<fCM) (3) 



where 4>^(r) is the envelope function of the n th state of 
a single carrier of type p, e is the electron charge, e is 
the dielectric constant of the QD material, and S n is the 
(pseudo) spin of state n. 

The elliptic disk shape of our model QD is symmetric 
under reflections about planes perpendicular to its main 
symmetry axes. Therefore, single carrier envelope func- 
tions are either odd or even under these reflections. The 
term s^rpn is even under the application of the same 
reflection for both r\ and r*2- Therefore, the parity of 



the integrand in Eq. ([3]) under such a 'double reflection' 
is determined by the parities of the envelope functions 
only. Whenever the integrand is odd under a 'double re- 
flection', the integral vanishes. We use these symmetry 
considerations in order to reduce the required computa- 
tion resources. 
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C. Electron- hole exchange interaction 

The zeroth order term in the long range EHEI , 
^jzMtjs U ©) ec l ua l s zero. Therefore, higher or- 

der terms in the long-range as well as the short-range 
exchange integral must be considere d 29 ' 33 . 
The pseudo-spin structure of the EHEI for the lowest en- 



ergy envelope functions is deduced from symmetry con- 
siderations (the method of invariants) 11 ' 20 . The SVEFA 
requires that the same considerations hold also for any 
other combination of envelope functions^ 9 -. Thus, we ex- 
press the electron- hole-exchange terms Cjffiij 3 i 4 as fol- 
lows: 



C: 



hehe 



a n,j2,nh,nj 3 ,nii* 






a; 



n J 2 ,nj 3 ,n?4 



^0 









a nj 2 ,nii,nj3,ni4 

a nj 2 ,Tlii ,nj 3 ,ni&* 
^2 






A.nj 2 ,ni 1 ,nj 3 ,ni 4 
^2 

\nj 2 ,ni 1: nj 3 ,ni 4 
^0 



(4) 



Where nik represents the index of the envelope 
function belonging to state number ik- The e-h 
pseudo spin base for the matrix are the functions: 

{i it), i w, i n), 1 w> }. 

The terms Ao and A2 are mainly affected by the short- 
range interactio n 20 ! 29 . This intra-unit-cell interaction is 
not sensitive to the details of the slowly varying envelope 
wavef unctions^ 9 -. Therefore, we assume that all the non- 
vanishing Aq 2 '™ 1 '™ 3 '™ 4 and A^ 2 '™ 1 '" 3 '™ 4 terms have the 
same values, Aq and A2, respectively. The values that 



we used were chosen such that the calculated X° spec- 
trum would fit the magneto-PL measured X° spectrum 12 . 
Since the short-range interaction is even under 'double 
reflections—, the symmetry considerations which aid in 
identifying the vanishing Coulomb integrals apply also in 
identifying the vanishing EHEI terms A and A 2 . 
The A™ 2 '™ 1 '™ 3 '™ 4 integrals are mainly affected by the sec- 
ond order terms in the expansion of the long-range inter- 
action, which are given b y 20 ' 29 



l A n a ,n 1 ,n„n 4 = j j ^ g^l^^^ ^^)^^)^^ (5) 



where n is a unit vector in the direction of f± — f 2 , 1 is the 
3x3 unit matrix, and is the valence-conduction band 
dipole matrix element. The dipole matrix element (/2) is 
related to the momentum matrix element (M) through 
the particle's mass and the energy difference between the 
dipole states 3 ^ (the bandgap energy E g ), 



m E q 



-M- 



(6) 



where the conduction-valence band momentum matrix 
elements are given by , 28 ' 31 



M T , t = = 



(7) 



M T4(l,t) = 2\/ m o^ P (l,(-)«,0) 



(8) 



where E p is the bulk material conduction-valence band 
interaction energ y 31 ' 35 and the spin quantization axis is 
chosen along the [001] (or z) direction. For compatibil- 
ity with the experimentally defined axes, we choose the 
major axis of the QD, believed to be along the [110] crys- 
tallographic axis^ as the x (or H) direction. 
Substituting Eqs. ([5115]) into Eq. © yields 



t n a ,ni,na,Tu _ 3e?h 2 E p f f ^my* \ (gi ~ ^) 2 - fri - x 2 ) 2 + 2i( Xl - x 2 ){yi - y 2 ) in , tfUnitf ^ 3 „ * 



2em Q E 2 g J J 9h { 2> (( Xl - x 2 )* + ( Vl - y 2 f + {z x - z 2 f)i ^ \ V9 e K D 

(9) 
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One can now, in principle, compute these integrals us- 
ing the single carrier envelope wavefunctions that were 
numerically obtained earlier. This approach demands 
a lot of computation resources in order to obtain reli- 
able accuracy. Therefore we choose to approximate the 
wavefunctions analytically, by using in-plane harmonic 
oscillator functions^ instead of the numerical ones (see 
Appendix |A| . With these approximations the 6 dimen- 
sional integrals are reduced into 5 analytical ones 3 -!. The 
non-analytical integral can be easily calculated numeri- 
cally. Alternatively, for a nearly round QD, this integral 
can be expanded into a power series in the aspect ratio of 
the model QD, from which only terms up to the linear or- 
der can be kept. This approach provides also important 
insight. The result of this derivation for A,' ' ' is 



A 



3y^e 2 h 2 E p (£ - 1) 



(10) 



where 1% is the characteristic length of the electron (Gaus- 
sian) wavefunction in the x direction (see Appendix [A"]l . 
/3 = 0.72 is the ratio between the characteristic length 
of the hole wavefunction to that of the electron, and 
£ = 0.96 is the length ratio between the short and long 
sides of the rectangle (the aspect ratio). 
For 1% = 72A, which gives the same s-p x energy separa- 
tion for the electrons as the numerical wavefunctions, we 
calculated the Ai terms that we list in Table ITTT In the 
table we also list the values that we could directly deduce 
from the measured fine structure splitting of the X° and 
the X -2 lines (A}' 1 ' 1 ' 1 and A}' 2 ' 1,2 , respectively). The 
agreement, as can be seen in the table, is remarkable. 
We note that the ratios a" 2 '™ 1 '™ 3 '"" 4 to A x ' ' ' can be 
quite large for small deviations from cylindrical symme- 
try. In particular, there are significant sign variations 
between the various terms. The expressions for these ra- 
tios as functions of (3 and £ (for |1 — £| <C 1) are also 
given in table HT1 



D. Optical transitions: Polarization selection rules 

The optical transition operator in the dipole approxi- 
mation is expressed as>^ 



(11) 



*,3 



Under the one-band SVEFA, the transition momentum 
vector pij is given by: 



(12) 



The momentum matrix elements Mg. ^ are given explic- 
itly by Eqs. (|7ll8l) . 

The rate of an optical transition^ centered at an energy 
e, for a certain polarization e is given by: 

Aane 



Te(e) 



ihrrvkc 2 



J2\ifW-P\i)\%, ei -e f Fi (13) 



TABLE II: Calculated, measured and estimated electron-hole 
exchange interaction terms (in fieV). The measured and esti- 
mated terms were used for calculating the PL spectra. The 
calculated ratios are given in terms of the hole/electron length 
ratio (3 and the aspect ratio £, for [1 — £| <C 1. 



Parameter 



A 
A 2 

A i.i.i,i 

A l,2,l,2 



A l,3,l,3 

A M,M 

A l,5,l,5 
A 2,l,2,l 
A 2,2,2,2 
A 2,3,2,3 
A 2,2,2,3 
A l,2.1,3 



Used in fit 



207 
21 
-25 
196 



-222 
-6.4 
232 
379 
119 
-12 
71i 
209i 



Calculated 



-15 
118 



-133 
-3.8 
139 
227 
71 
-7.4 
42i 
125i 



Ratio to A 1 



1,1,1,1 



2g-l 



a A 5-2 

7JI 



1+/3 2 €-1 
f3 4 41-6g 
(l+(3 2 ) 2 16 

p 2 1 r ff 4 61g-45 

T+^e-1 ^ (1+/3 2 ) 2 32(C-1) 
1 2g-l 

T+W 1 €-i 

li 2 61|-45 
(1+/3 2 ) 2 16(4-1) 
/3 2 41-6g 
(1+/3 2 ) 2 16 
■ /3 2 . 9-g 

l (i+0 2 y i i6(e-i) 



a 2 



l+i 



'l+W 1 2(5-1) 



where a = w is the fine structure constant, and n 
is the refraction index of the QD material. The indices i 
and / run over all initial states \i) and final states |/). Ej 
(ef) is the energy of the initial state \i) (final state |/)). 
Fi is the population probability of the initial state 
For the bright neutral exciton transitions we calculate 



\m-p\x° H ) 



= \(0\y-P\X° v )\ 2 = 



m a E p 



1.44 (14) 



The other two transitions from the bright states and the 
transitions from the 'dark' excitonic states, vanish. As- 
suming equal population probabilities for the bright and 
dark X° states, we get a total X° rate of (0.78ns) _1 , in 
agreement with the measured lifetime 7 . The calculated 
rates of all other optical transitions are given in units of 
this total X° rate. 

For example, in Fig. U we present a diagram of the calcu- 
lated many carriers energy levels and the optical transi- 
tions between them, which lead to the spectrum resulted 
from excitonic transitions in a triply negatively charged 
QD(X- 3 ). 

As can be deduced from Fig.dJ the X~ 3 line results from 
three initial levels (each doubly Kramres degenerate). 
These levels contain mainly the following single carrier 
states: one s-shell hole, two s-shell electrons one p x and 
one p y electron, where the p-shell electrons are in their 
triplet configurations. These open shells configurations 
are the lowest in energy, since the energy difference be- 
tween the p x and p y single electron states is smaller than 
their exchange interaction. The degeneracy between the 
triplet configurations is removed by the EHEI with the 
hole. Our experimental data can only be explained with 
these open p-shells occupation in mind (see below). 
The final states are mainly composed of three single elec- 
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tron states: one in the s-shell, one in the p x shell and one 
in the p y shell. The expected eight fold degeneracy is 
partially removed by the electron - electron exchange in- 
teraction, which leave a four fold degenerate ground state 
(we found no experimental evidence for an anisotropic e-e 
exchange interaction 9 which would have further reduced 
this degeneracy). The calculated optically allowed tran- 
sitions between the initial and the final states and their 
polarization selection rules are given in Fig.[4j The high- 
est energy transition is finely structured from three lines 
with intensity ratios of approximately 3:2:1. These in- 
tensities were previously deduced using a simple model, 
by Urbaszeck et a& 

If the ground p-shell was occupied by two electrons, the 
X~ 3 transitions would have generate a single almost un- 
polarized spectral line, very similar to that due to the 
X -1 transitions. This is in clear contradiction with the 
measurements presented in Fig. [3l Fig. [8] and previous 
measurements on similar QDs'i'21. 

Another example is provided in Fig. [5] where we show 
the levels' diagram and optical transitions, which result 
in the biexciton recombination in a doubly negatively 
charged QD (XX~ 2 ). 

Here the initial states are mainly composed of the same 
single carrier states as the initial states of the X -3 , ex- 
cept for the addition of one more s-shell hole. Unlike the 
X" 3 , the paired s-shell holes do not remove the degen- 
eracy of the triplet configuration of the p-shcll electrons. 
Similarly, the final states are mainly composed of the 
same single electron states as the final states of the X~ 3 , 
except for an additional one s-shell hole. The EHEI be- 
tween the unpaired hole and the electrons, completely 
removes the degeneracy between the electron states. As 
a result, there are eight low energy states (the energy 
differences between the lowest and between the highest 
pairs of states are too small to be noted) to six of which 
optical transitions are allowed. Similar to the case of the 
X -3 , the optical spectrum is finely structured from three 
pairs of lines with total intensity ratios of approximately 
3:2:1 as previously deduced by the simple considerations 
of Urbaszeck et al^. A major difference between the two 
examples is in their polarization selection rules. In the 
first case (X -3 ) the total spin is half integer and Kramers 
degeneracy prevails. Therefore only partial linear polar- 
ization is expected. In the second case (XX -2 ) the total 
spin is an integer. In this case, full linear polarization is 
expected, just like in the case of the neutral single ex- 
citon. Indeed, the intermediate pair of spectral lines is 
fully polarized along the QD's primary axes. 
Similar diagrams for less complicated transitions were 
discussed in previous work a 11 ' 21 ! 39 . In these works sim- 
pler models were used. These models are sufficient only 
when the EHEIs are much smaller than any other inter- 
action. 

In Figs. EJa) and [6jc) we show the calculated spectrum 
for various charge states. The corresponding H-V po- 
larization projections are shown in Figs. [6]Jb) and[6jd), 
respectively. 
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FIG. 4: Schematic description of the calculated many carri- 
ers energy levels, and their spin wavefunctions, which optical 
transition between them result in the X -3 spectral lines. Each 
transition rate (in units of the total X° rate) and its degree 
of polarization are indicated (if absent the transition is un- 
polarized). Only one of the two Kramers states is shown for 
each level (for notation see Appendix[B}. The number before 
each component indicates its amplitude. Components which 
are irrelevant to the polarization degree and have amplitudes 
below 0.1 are not included. 



Within our simple, one band model, the calculated D-D 
(and of course the R-L) projections vanish, and therefore 
they are not shown. For the calculations, equal proba- 
bilities for excitons and biexcitons in all charge states 
were assumed^ -. In the calculations, only initial configu- 
rations within 1 meV (compatible with the experiments' 
temperature) above the ground state were considered. 
The calculated lines are convoluted with 50^eV broad 
Gaussians, to account for the spectral diffusion^. In the 
calculation of the polarization projections, a constant 
background of 3.5% of the maximal intensity is added 
to both cross-linearly polarized spectra. This is done 
in order to mimic the effect of background noise on the 
measured spectra (see sect. HT)) . 



IV. COMPARISON BETWEEN 
EXPERIMENTAL MEASUREMENTS AND 
MODEL CALCULATIONS 

In Fig. [7] we compare between the measured and cal- 
culated spectral positions of various lines. 
We note that the spectrum produced by our simple model 
correlates with the experimentally measured one in the 
energy order of the various spectral lines. The calculated 
energy differences between the various lines, however, are 
in most cases smaller than those measured. Specifically, 
the calculated energy differences between the exciton and 
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FIG. 5: Schematic description of the calculated many carri- 
ers energy levels, and their spin wavefunctions, which optical 
transition between them result in the XX -2 spectral lines. 
Each transition rate (in units of the total X° rate) and its 
degree of polarization are indicated (if absent the transition 
is unpolarized) . All states are shown (for notation see Ap- 
pendix [Bj. The number before each component indicates its 
amplitude. Components which have amplitudes below 0.1 are 
not included. 



the biexciton and between the positive and negative tri- 
ons (X +1 and X -1 , respectively) are smaller than the 
measured ones. This is probably a consequence of the rel- 
ative simplicity of our single band model^ and the lack 
of information about the exact shape strain and com- 
position of the QDs. With our model's limitations we 
found it hard to simultaneously fit the biexciton binding 
energy and the difference between the positive and neg- 
ative trion transitions. 

In Fig. [5] we compare the measured and calculated po- 
larized fine structure of various spectral lines, while In 
Fig. [9] we compare the measured and calculated linear 
polarization spectra for these spectral lines. 
We note in Figs.[5]and[5]that the measured fine structures 
are reproduced quite nicely by our model calculations. In 
particular, the calculated number of fine structure com- 
ponents, their relative intensities and their polarizations 
correlate with the measured values. 
On the other hand, while the calculated polarization 
spectra are always polarized along the H-V axis of the 
poincare sphere, the measured ones are sometimes ro- 
tated. Few specific lines (see below) are polarized along 
the (V+D)-(H+D) axis. 

The calculated energy differences between the fine struc- 
ture components of a given spectral line are sometimes 
larger than the measured values. Particularly, the cal- 
culated fine structure splittings between the components 
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FIG. 6: Calculated PL spectra (a), (c), and their polarization 
projections on the H-V axis (b), (d), for various single QD 
excitonic transitions. Vertical dash lines at various spectral 
lines are drawn to guide the eye. 



of the X -3 line, and the calculated splitting between the 
unpolarized and polarized doublets of the X -2 line are 
larger than the measured ones. We believe that this 
maybe a consequence of the An dependencies on the en- 
velope wave functions, which are neglected in our model. 
In the absence of polarization memory (which requires 
quasi-resonant polarized excitation ») , and in the absence 
of magnetic field, the theory predicts that the spectral 
lines can only present linear polarizations. This is what 
we observed experimentally as well. In the theoretical 
model, the linear polarization can only be oriented along 
the main axes of the QD, which are usually along the 
crystalline directions [110] and [110]^, which we denote 
by H and V respectively. In the experiment, however, 
we found that some spectral lines are polarized along 
other directions. The polarization appears in three sets 
of orthogonal axes: The measured polarization of the 
neutral exciton line is indeed along the H-V axis of the 
Poincare sphere. A few other lines, notably the neutral 
biexciton, the doubly charged biexciton and the triply 
charged exciton, are also polarized along this axis. Their 
degree of polarization is somewhat smaller than that of 
the X° lines due to the unpolarized spectral background 
that they ride on. Few other spectral lines are polarized 
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FIG. 7: A comparison between measured and calculated spec- 
tral positions of a few lines. The dashed line is the equality 
line. The size of the markers represent the experimental error. 



along an axis which is rotated clockwise by 135 degrees 
relative to the H-V axis of the Poincare sphere. 
This polarization axis, which roughly coincides with the 
-^(V+D) ([120]) and the ^(H+D) ([210]) crystalline 
directions, appears only in lines associated with config- 
urations which contain one unpaired p x carrier (either 
electron or hole): X" 2 , XX" 1 and XX+ 1 . Careful in- 
spection of the unevenly polarized spectra of the X -2 
line presented by Ediger et al21 leads to the same con- 
clusion. Unfortunately, they did not fully measured the 
polarization state of the line. 

For other configurations, which contain only s carriers 
or either closed shells or two unpaired p carriers (p x and 
p y ), the polarization is along the H-V axis. In addition, 
we sometimes observe lines which are polarized along the 
D-D axis as well. Such a spectral line is seen in Fig.^d), 
at energy of -2.6 meV relative to the X° line. This rela- 
tively weak doublet, may result from pair recombination 
in doubly negatively charged QD as deduced from its 
voltage dependence. 

These novel observations are not reproduced by the sim- 
ple model that we present, probably since the model ig- 
nores the underlying crystal. For example, in the single 
band model the polarization due to recombination of a 
p x shell e-h pair is the same as that of a p y pair, which 
is clearly not the case. 

We rule out another possible explanation 42 expected in 
a nearly p-shell degenerated QD. In such a QD, where 
the p x -Py splitting is comparable to or smaller than the 
EHEI, configurations containing an unpaired p x carrier 
and those containing an unpaired p y carrier are mixed 4 ^. 
This mixing can indeed lead to recombination in linear 
polarizations along axes different than the primary axes 
of the QD. This is not the case here, since for nearly de- 
generate QDs, the fine structures of the X -2 , XX -1 and 
XX +1 transitions should contain twice the number of fine 
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FIG. 8: (color online) Measured (top panel in each pair) and 
calculated (bottom panels) high resolution polarization sensi- 
tive PL spectra of various spectral lines. The solid (dashed) 
black line represents H (V) polarized spectrum while the solid 
(dashed) pink line represents the V+D (H+D) polarized spec- 
trum. 



structure components than that actually observed. 
One straight forward way to include the lattice will be 
to solve the single electron problem using a multi-band 
approac h 31 ' 43 . 

Such a model is absolutely necessary for calculating the 
polarization selection rules in highly positively charged 
QDs. There, the complicated nature of the p-shell holes 
is not likely to be captured by a one band model. The 
polarization degree of the highly positively charged QD 
emission lines (identified by their voltage dependence) 
that we measured was marginal, while our single band 
model yields polarizations similar to those of the nega- 
tively charged QD lines. We believe that this discrepancy 
results from the inadequacy of the single band model. In 
addition, we note that the highly positively charged QD 
lines were all measured under large electrostatic fields, 
which for now, were not considered in our model. 
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are favorably compared with the measured polarization 
sensitive PL spectra. 

However, while our model can only reproduce polariza- 
tions along the main axes of the QD, the experimental 
data display polarizations oriented along other direc- 
tions as well. This indicates, probably, that the one 
band based model is too simple to describe this novel 
observation. 
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APPENDIX A: THE ANALYTICAL 
WAVEFUNCTIONS USED IN THE 
CALCULATION OF THE Ai INTEGRALS 

As model functions we used the following functions: 
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FIG. 9: (color online) Measured (top panels) and calculated 
(bottom panels) linear polarization spectra of various spectral 
lines. Black (pink) solid lines represent projections along the 
H-V ((V+D)-(H+D)) axis of the poincare sphere. 



V. SUMMARY 

We presented detailed polarization sensitive spec- 
troscopy of single QDs in various charge states. Wc 
developed a many-carrier model based on single band, 
envelope function approximation, which includes the 
isotropic and anisotropic electron-hole exchange in- 
teractions for the analysis of the measured data. We 
calculated the PL spectrum with its fine structure and 
polarizations for the exciton and biexciton optical tran- 
sitions in neutral, singly positively and singly, doubly 
and triply negatively charged QDs. The calculations 



The kets of the form \n x , n y ) p stand for the 2D elliptic 
harmonic oscillator functions: 
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where p is the charge carrier index (either 'e' or '/i'), 
lp ^ is a characteristic length along the x(y) direction, 
n x (y) is t ne quantum number associated with the x(y) 
direction, and H n ^ y) is the Hermite polynomial of order 
n x (y)- The aspect ratio £ and the hole/electron length 

ratio (3 are defined as: £=jf=pr j P — $ = ^t- 



APPENDIX B: THE NOTATION USED FOR THE 
STATE VECTORS 

In table Mil we list the symbols and abbreviations used 
in writing the spin state vectors. 
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